AD-A220  735 


A  THEORETICAL  MODEL  FOR 
AIRBORNE  RADARS 

by 

D.  Faubert 


DEFENCE  RESEARCH  ESTABLISHMENT  OTTAWA 

REPORT  NO.  1017 


Canada 

DiSTr;.:-'':-.-  ■  ■  '  ■;  :  a 

A;-'.  ■  ■'  •  .l.-r.se; 

...  ^  ^  -iicd 

November  1989 
Ottawa 

90  (!■' 

National  Defense 
Defence  nadonale 


A  THEORETICAL  MODEL  FOR 
AIRBORNE  RADARS 


D.  Faubert 

Airborne  Radar  Section 
Radar  Division 


DEFENCE  RESEARCH  ESTABLISHMENT  OTTAWA 

REPORT  NO.  1017 


PCN 
21  LAI  2 


November  1989 
Ottawa 


ABSTRACT 


This  work  describes  a  theoretical  model  for  the  simulation  of  airborne  (or 
spaceborne)  radars.  This  model  can  be  used  to  simulate  many  types  of  systems 
including  Airborne  Intercept  and  Airborne  Early  Warning  radars,  Airborne  Missile 
Approach  Warning  Systems  etc.  We  derive  a  theoretical  expression  for  the  average 
Signal— to— Noise  ratio  at  the  output  of  the  signal  processor.  This  can  be  used  to  obtain 
the  average  performance  of  the  radar  without  employing  Monte  Carlo  techniques.  The 
model  makes  provision  for  a  waveform  without  frequency  modulation  and  for  one  with 
linear  frequency  modulation.  The  waveform  may  also  use  frequency  hopping  for 
Electronic  Counter  Counter  Measures  or  for  clutter  suppression.  The  model  can 
accommodate  any  type  of  encounter  including  air— to— air,  air-to-ground  (look— down) 
and  rear  attacks.  It  can  simulate  systems  with  multiple  phase  centers  on  receive  for 
studying  advanced  clutter  or  jamming  interference  suppression  techniques.  To 
demonstrate  the  validity  and  the  capability  of  the  model,  a  baseline  Airborne  Intercept 
radar  system  is  simulated. 


RESUME 

Ce  travail  decrit  un  modele  theorique  pour  la  simulation  des  radars  aeroportes 
(ou  deployes  dans  I’espace).  On  pent  I’utiliser  pour  simuler  plusieurs  genres  de  systemes 
incluant  les  radars  d’interception  aerienne,  les  radars  de  veille  a  longue  portee,  les 
systemes  de  surveillance  d’attaque  par  missiles  etc.  Le  modele  calcule  le  rapport  signal 
sur  bruit  moyen  du  systeme  apres  traitement  du  signal.  II  peut  simuler  une  onde 
porteuse  sans  modulation  de  frequence  ou  avec  modulation  lineaire  de  la  frequence. 
L’onde  porteuse  peut  changer  de  frequence  period! quement  pour  tirer  avantage  des 
techniques  de  contre— centre  mesures  electroniques  ou  de  reduction  du  clutter.  Le 
modele  s’applique  a  des  engagements  air— air,  air— sol  et  d’attaque  par  I’arriere.  II  peut 
decrire  des  recepteurs  ayant  plusieurs  centres  de  phase.  Ce  genre  de  recepteur  est 
quelques  fois  utilise  pour  reduire  I’lnterference  cause  par  le  clutter  ou  le  brouillage 
electronique.  On  presente  quelques  resultats  ayant  trait  a  un  radar  d’interception 
aerienne  pour  demontrer  I’utilite  et  la  validite  du  modele. 

iii 


EXECUTIVE  SUMMARY 


The  requirement  for  evaluating  the  performance  of  airborne  radars  often  arises  in 
the  Department  of  National  Defense.  This  report  is  the  outcome  of  an  effort  to  develop 
a  general  purpose  model  for  the  simulation  of  these  systems.  It  can  be  used  to  study 
multi— mode  Airborne  Intercept  radars  (AI)  deployed  on— board  high  performance  fighter 
aircraft,  Missile  Approach  Warning  Systems  deployed  on  helicopters,  propeller  and  jet 
aircraft.  Airborne  Early  Warning  radars  deployed  on  surveillance  aircraft  like  the 
AWACS,  etc.  The  model  also  applies  for  Space  Based  Radars  (SBRs). 

In  the  theory,  we  derive  expressions  for  the  signals  found  in  the  radar  receiver. 
Hence  with  minor  modifications,  the  model  could  be  used  to  generate  signals  for  input  to 
a  real  radar  signal  processor. 

To  ensure  the  widest  possible  applicability,  the  model  accounts  for  various 
system  configurations.  It  has  provision  for  various  types  of  waveforms.  The  radar 
carrier  frequency  can  also  change  periodically  to  counter  jamming  or  to  reduce  clutter 
interference.  The  model  also  applies  to  radars  equipped  with  more  than  one  antenna,  as 
radars  of  this  type  perform  better  in  the  presence  of  jamming  and  clutter.  Finally,  the 
model  geometry  accommodates  air— to— air,  air-to-ground  and  rear— attack  encounters. 

The  first  part  of  the  report  contains  the  equations  for  the  signal  received  from  a 
point  target.  The  signal  from  a  distributed  target,  like  ground  clutter,  is  an  extension  of 
the  expressions  for  the  point  target.  The  calculations  also  deal  with  the  contribution  of 
the  receiver  thermal  noise.  Finally,  an  expression  for  the  average  value  of  the  combined 
interference  of  clutter  and  noise  is  derived.  This  allows  the  evaluation  of  the  average 
performance  of  the  radar  without  using  time-consuming  Monte  Carlo  techniques. 


V 


A  baseline  radar  system  serves  as  a  test  bed  to  demonstrate  the  validity  and  the 
capability  of  the  theory.  The  model  highlights  the  impact  of  clutter  and  white  thermal 
noise  on  target  detection. 

The  model  has  a  wide  range  of  applications  for  DND.  It  could  be  used  to 
evaluate  Airborne  Intercept  and  Airborne  Early  Warning  radars,  Missile  Approach 
Warning  Systems  etc.  It  is  also  of  interest  to  the  TTCP  Panel  KTP— 5  which  is 
concerned  with  the  detection  of  low  radar  cross  section  targets. 


VI 
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1.0  INTRODUCTION 


Modern  military  aircraft  make  extensive  use  of  radar  systems  for  target 
acquisition.  As  an  example,  the  Canadian  CF— 18  aircraft  is  equipped  with  the  Hughes 
AN/APG— 65  coherent  radar.  It  uses  pulse  Doppler  techniques  to  counter  clutter  and 
jamming.  This  multi— mode  Airborne  Intercept  (AI)  radar  can  operate  in  a  number  of 
ways  to  deal  with  a  wide  variety  of  targets  under  various  operational  conditions.  The 
importance  of  radars  in  many  airborne  weapon  systems  used  by  the  Canadian  Forces 
dictates  a  detailed  understanding  of  their  operation  and  capabilities. 

In  support  of  the  Canadian  Forces  (CF),  the  Radar  Division  of  DREO  is  studying 
airborne  radar  systems  under  Tasks  DFTEM— 102.  This  report  describes  the  first  phase 
of  this  work  consisting  in  modeling  the  operation  of  airborne  radars. 

Section  2.0  derives  the  airborne  radar  model.  Section  3.0  describes  the  computer 
simulation  implementing  the  theoretical  model.  This  section  also  defines  the  baseline 
airborne  radar  used  to  demonstrate  the  capability  of  the  simulation.  Section  4.0 
summarizes  the  findings  of  this  report. 

2.0  THEORETICAL  MODEL 

The  theory  developed  for  modeling  airborne  radars  (ABR)  is  an  extension  of  the 
one  derived  in  Refs,  l-^  for  Space  Based  Radars  (SBR).  It  is  a  so-called  video 
simulation  which  models  the  signals  found  in  the  video  section  of  the  radar.  The  results 
of  Refs.  1-4  do  not  apply  directly  to  the  airborne  radar  problem.  For  SBRs,  it  is 
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appropriate  to  neglect  second  order  phase  terms  but,  as  will  be  seen  later,  they  may  be 
significant  for  ABRs.  The  SBR  model  did  not  simulate  pulse  compression  waveforms. 
Eclipsing  does  not  normally  constitute  a  problem  for  SBRs  but  occurs  for  ABRs. 
Finally,  the  present  work  allows  for  a  frequency  stepped  waveform  and  the  use  of  guard 
antennas,  features  not  included  in  the  SBR  model. 

Table  I  describes  the  capabilities  of  the  simulation.  Although  some  of  these 
options  are  not  available  in  the  equipment  used  by  the  CF,  they  are  included  in  order  to 
allow  for  the  investigation  of  more  advanced  systems. 


TABLE  I 

AIRBORNE  RADAR  SIMULATOR  OPTIONS 


CONFIGURATIONS  Pulse  Doppler  Radar 

Displaced  Phase  Center  Antenna 
Multiple  Phase  Centers 
Guard  Antennas 

WAVEFORMS  No  Frequency  Modulation 

Linear  Frequency  Modulation 
Frequency  Hopping 
Clutter— lock  circuitry 

MODES  OF  OPERATION  Air-to-Air 

Air— to— Ground  (Look— Down) 
Rear  Attacks 


Table  I  shows  the  various  Airborne  Radar  Simulator  options.  For  the  four 
configurations  listed  in  Table  I,  the  location  of  the  radar  transmit  phase  center  is  in  the 
middle  of  the  antenna.  In  the  conventional  Pulse  Doppler  Radar  (PDR)  configuration, 
the  locations  of  the  receive  and  the  transmit  phase  centers  coincide.  In  the  Displaced 


2 


Phase  Center  Antenna  (DPCA)  configuration,  there  are  typically  two  or  three  receive 
phase  centers.  Their  locations  are  such  that  the  motion  of  the  airborne  platform  does 
not  degrade  clutter  cancellation.  In  the  Multiple  Phase  Centers  configuration,  the  radar 
uses  a  large  number  of  receive  phase  centers.  Appropriate  weighting  of  the  received 
signals  by  the  antenna  beamformer  reduces  jamming  interference  and  possibly  clutter. 
In  the  Guard  Antenna  configuration,  the  signal  on  the  main  antenna  is  weighted  with 
the  signal  from  one  or  more  separate  guard  antennas.  This  achieves  jamming  and 
clutter  reduction  through  antenna  sidelobe  nulling  or  blanking. 

The  waveform  is  a  burst  of  radio  frequency  pulses  without  frequency  modulation 
(NFM)  or  with  linear  frequency  modulation  (LFM)  for  pulse  compression.  For  the 
detection  of  fast  moving  targets,  the  usual  choice  is  the  NFM  waveform  with  a  high 
pulse  repetition  frequency  (HPRF).  The  LFM  mode  normally  goes  with  a  low  pulse 
repetition  frequency  (LPRF)  for  range  unambiguous  operation  and  with  a  medium  pulse 
repetition  frequency  (MPRF)  for  the  detection  of  ground  targets  (look-down  mode). 
The  carrier  frequency  of  the  waveform  can  change  from  pulse  to  pulse  to  improve  the 
performance  of  the  radar  against  clutter  and  jamming. 

In  the  following  sections,  a  model  for  the  simulation  of  ABRs  is  developed.  This 
model  includes  computation  of  the  signals  for  a  point  target,  clutter  and  the  receiver 
thermal  noise.  Finally,  the  average  Signal-to-Noise  ratio  at  the  output  of  the  radar 
signal  processor  is  derived.  The  model  allows  one  to  study  the  average  performance  of 
airborne  radars  without  having  to  use  time-consuming  Monte  Carlo  techniques. 
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2.1  The  Geometry 


The  calculation  of  the  radar  signals  employs  a  geocentric  frame  of  reference  shown 

in  Fig.  1.  The  earth  is  assumed  to  be  circular  with  a  radius  equal  to  The  value  of  r^ 

normally  equals  6378  km,  but  can  be  multiplied  by  4/3  or  any  other  appropriate  factor 

to  account  for  the  effect  of  the  refraction  of  the  atmosphere.  In  Fig.  1,  T  ,  R-  and  C 

represent  respectively  the  locations  of  the  transmit  phase  center,  the  i*'^  receive  phase 

center  and  a  point  scatterer  (target  or  clutter).  The  distance  between  T  and  R-  equals 

d/2.  Let  r  [t]  and  r  [t]  be  the  vectors  describing  the  position  of  T  and  C,  and  i  be  a 
SC  s 

unit  vector  along  the  line  joining  T  and  R..  It  should  be  noted  that  throughout  this 

report,  square  brackets  always  indicate  the  functional  dependence  of  a  variable  on  the 

parameter  within  the  brackets.  The  radar  (both  T  and  R.)  and  the  target  are  moving 

at  velocities  v  [t]  and  v  [t]  respectively.  The  velocity  vector  v  [t]  includes  the 

contribution  from  the  airspeed,  the  yaw,  the  roll  and  the  pitch  of  the  airplane.  The 

velocity  of  the  radar  is  not  necessarily  equal  to  the  velocity  of  the  center  of  mass  of  the 

aircraft.  This  is  immaterial  in  our  calculations  since  we  model  the  radar  independently 

from  the  aircraft.  The  time  elapsed  between  the  transmission  and  the  reception  of  the 

waveform  is  of  the  order  of  lO'^s  for  a  typical  ABR.  Hence,  v  [t]  and  v  [t]  can  be 

s  c 

considered  constant  during  the  measurement  and  equal  to  their  values  at  t=0. 

Without  loss  of  generality,  one  can  specify  that  at  t=0,  T  rests  on  the  Z  axis  of 

b 

the  calculation  reference  frame.  At  this  time,  a  burst  of  rectangular  RF  pulses  is 
transmitted  from  the  radar  transmitter.  The  field  sensed  by  is  a  collection  of  delayed 
replicas  of  the  transmitted  waveform,  backscattered  by  the  targets  within  the  field  of 
view  of  the  radar.  The  delay  between  transmission  and  reception  of  the  waveform  is  a 
non-linear  function  of  time  which  is  the  subject  of  the  next  sub— section. 
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propagation  delay,  A.,  is  a  function  of  tg,  and  tj  and  is  given  by 


-  (<■(.  tg)  +  (tj  t^)  -  t.  tg  (1) 

The  term  (t^— tg)  can  be  evaluated  by  using  the  following  equation; 

where  c  is  the  speed  of  light  and  where  j  |  stands  for  the  norm  of  the  vector.  The 
vectors  r^  [t^]  and  r^  [tg]  in  Eq.  2  can  be  evaluated  from  the  following  relations; 


=  'cl^l  + 

(3) 

and 

'si'ol  =  'sl°l  +  ‘o^sl°l 

(4) 

In  order  to  simplify  the  notation,  we  omit  the  functional  dependence  on  time  of  the 

various  parameters  involved  in  the  calculations  when  t=0. 

With  this  convention  and 

the  help  of  Eqs.  3  and  4,  Eq.  2  becomes: 

=  lr  +tv  — r  -  t^v  1 
'cO''  'c  cc  s  Os' 

(5) 

-♦ 

Let  R(t)  be  the  vector  joining  T  to  C; 

o 

R(t)  =  r  Jt]  -  ? Jt] 

(6) 
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Expanding  Eq.  5  and  using  Eq.  6,  one  obtains 


(l<.-‘o)  =  {|R|^  +  +  l‘cVVsl^r^^  (’’) 

To  further  simplify  the  notation,  we  represent  the  norm  of  a  vector  by  omitting 
the  arrow  (  )  on  top  of  the  corresponding  vector. 


Then  we  expand  Eq.  7  into  a  Taylor  power  series 


=  R  +  VR-v^-tf^R-v^+ 


c  c  Os 


2 


c  c  0  s 


- 7? 


2R 
2  2 


+  ... 


(8) 


For  SBR’s,  the  quadratic  and  higher  order  terms  of  Eq.  8  (in  italics)  are  small 
compared  to  the  wavelength  of  the  radar.  This  is  because  for  SBRs,  the  range  to  the 
closest  target,  R,  is  hundreds  or  even  thousands  of  kilometers  and  is  typically  an 
order  of  magnitude  longer  than  for  ABRs.  By  using  the  worst  case  values  given  in  Table 
II,  it  is  readily  apparent  that  for  ABRs,  the  second  order  terms  of  Eq.  8  can  be 
significant.  Neglecting  cubic  and  higher  order  terms  which  are  negligible,  Eq.  8 
becomes: 


c(tc-to)  =  R  +  t^R'V^  -  t^R-v^+ 


‘c''c  -‘o''sl 


0' 


1 

m 


tcR-v 


2R 


(9) 
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TABLE  n 


WORST  CASE  VALUES  FOR  ABR  CALCULATIONS 


Maximum  platform  velocity  (v  ) 
Maximum  target  velocity  (v^) 
Maximum  radial  velocity  (V^) 
Maximum  normal  velocity  (V^) 
Target  maximum  range 
Target  minimum  range 
Maximum  time  elapsed  (tg,  t^,  t|) 


Mach  1  (approximately  300  m/s) 
Mach  1 

Mach  2  (approximately  600  m/s) 
Mach  1 


100  km 

100  m 

_o 
10  ^s 


Eq.  9  can  be  simplified  further.  The  value  of  (t(«~tQ)  is  approximately  equal  to 
2R/c.  Hence  for  small  values  of  R  ,  (t^-tg)  is  small  and  it  can  be  assume  that,  for  the 
quadratic  terms  only,  For  large  values  of  R,  the  approximation  breaks  down  but 

in  this  case,  the  value  of  the  quadratic  terms  is  small  due  to  the  R  ^  coefficient.  The 

_3 

error  made  by  using  this  approximation  is  of  the  order  of  10  A^.  Hence  Eq.  9  becomes: 


2R 


(10) 


The  relative  velocity  between  the  radar  and  the  target. 


V,  is 


(11) 
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and  the  magnitude  of  the  radial  relative  velocity,  V^.,  is 

Vj.  =  R.V  (12) 

The  magnitude  of  the  normal  component  of  the  relative  velocity  between  the  radar  and 
the  target,  V^,  is: 


(13) 


Using  Eq.  11—13,  the  propagation  delay  from  the  radar  to  the  target  can  therefore  be 
written  in  a  more  compact  form  as 


c(tc-to)  -  R  +  t^R-v^  -  t^R.  Vg  +  ^  (14) 

In  Eq.  14,  the  radial  component  of  the  relative  velocity  between  the  radar  and  the 
target  results  in  a  linear  phase  shift  of  the  signal  (Doppler  effect).  On  the  other  hand, 
the  normal  component  produces  a  quadratic  phase  term  (Ref.  5). 

The  next  propagation  delay  which  must  be  evaluated  in  Eq.  1  is  c(tj-t^).  This 
term  can  be  calculated  by  use  of  the  following  relationship: 

c(tj-t^)  =  I  rj[tj]  -  ? Jt  J I  (15) 
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The  vector  r-(t)  can  be  written: 


rjW  =  rjt]  +  {d/2)i  +  t-v^  +  (d/2)i  (16) 

Following  the  procedure  used  to  derive  Eq.  14  and  dropping  third  and  higher  order 
terms,  one  obtains  after  some  tedious  but  straightforward  manipulations; 


c(tj-t^)=R  -  tjR-Vs  +  t^R-v^-^R-i - V- i  + 


dt :  .  t:  „  ,2 

-V^  +—  (17) 

2R  2R  ^  8R 


The  total  time  elapsed  from  transmission  to  reception  can  be  found  by  replacing 
Eqs.  14  and  17  in  Eq.  1.  However,  in  order  to  simplify  the  calculations  of  the  radar 
video  signal,  it  is  preferable  to  express  Aj  in  terms  of  a  single  time  variable,  t-  Hence, 
from  Eq.  14,  we  can  express  ctg  as; 


ctQ  - 


2  2 


R-v. 


1- 


-1 


(18) 


Expanding  Eq.  18  in  a  Taylor  power  series  and  neglecting  again  third  and  higher  order 
terms,  it  becomes 


R  +  -  c)  + 


c 


+ 


(19) 
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Using  Eq.  17,  we  can  compute  as  a  function  of  t-: 


(R.i)+a - t 

8R  ^ 


c- Vj+ 


(20) 


Using  Eqs.  19  and  20,  one  can  finally  compute  A.  from  Eq.  1.  Neglecting  again 
the  third  and  higher  order  terms,  the  result  for  the  overall  propagation  delay  is 


r  .R*  1  .  d  , 

~  +  ME  + 

c 


2V-(d/2R)V‘  i 


+ 


uftf 

n  i 
~1E~ 


(21) 


Eq.  21  is  identical  to  the  equations  obtained  in  Refs.  1,  4,  except  for  the  terms  in 
italics  which  are  second  order  effects  neglected  for  SBRs. 


Although  the  calculations  leading  to  Eq.  21  are  rather  tedious,  the  equation 
conveys  a  clear  physical  meaning.  The  term  2R/c  is  the  gross  delay  time  required  for 
the  RF  pulses  to  travel  back  and  forth  to  the  target.  The  term  2RV^/c2  accounts  for  the 
delay  resulting  from  the  distance  travelled  by  the  target  during  the  gross  delay  2R/c. 
The  term  -{d/2c)R-i  is  the  time  required  to  reach  the  i  phase  center.  This  can  be 
shown  by  referring  to  Fig.  2,  where  ARj  is  the  distance  between  points  A  and  R^; 


SR.  =  (d/2)  cos[^] 

=  —(d/2)  cos[(x/2)-|-tf] 

=  -(d/2)R-i  (22) 
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Therefore  — (d/2c)R*i  is  the  time  to  reach  the  i*^”  phase  centre.  The  term  d'^/8Rc  is  a 

2 

rather  small  second  order  effect  (d  /8R  ~  4>«10'2  Ac).  The  term  2tjV^  is  the  Doppler 

-♦ 

shift  experienced  by  the  RF  radiation  received  by  R-.  The  term  — tj(d/2R)V- i /c  is  a 
second  order  effect  arising  because  the  radial  velocity  of  the  target  and  hence  its  Doppler 

shift  is  a  function  of  the  location  of  R-.  This  can  be  explained  by  referring  to  Fig.  2. 

-*1  *♦ 

Let  Vj.  be  the  radial  velocity  of  point  C  relative  to  R-  and  Rj^  be  a  line  joining  Rj  and  C. 

We  thus  have: 

Vl  ~  (v  — V  )-R. 
r  ^  c  S’'  1C 

nc 

=  Vj.-(d/2R)  v.i  (23) 

2  2 

Finally,  the  term  V^tj/Rc  is  a  second  order  phase  term  caused  by  the  component  of  the 

■*  R 

motion  of  the  target  which  is  normal  to  V. 


Toward  C 

Fig.  2  —  Location  of  the  phase  centers  and  the  target 
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To  facilitate  the  derivation  and  the  interpretation  of  the  expressions  for  the  radar 
signals,  Aj[tj]  is  rewritten  by  grouping  together  the  various  powers  of  tj  and  dropping 
the  subscript  i  of  t-: 


where 


A.[t]  =  <5o  + 


,  _2R  ^^-^r  jR*i  ,  d^ 


^1  = 


2V^-(d/2R)V-i 


v2 

^o  =  — 

^  Rc 


(24) 

(25) 

(26) 

(27) 


The  radar  signal  received  from  a  target  is  a  replica  of  the  transmitted  signal,  delayed  by 
the  propagation  time.  The  modeling  of  the  transmitted  signal  is  the  subject  of  the 
following  section. 


2.3  The  Transmitted  Radar  Signal 


The  radar  transmits  a  waveform  described  by 

u[t]  =  gj  exp[juJ^t]  •  i/[t]  (28) 


where  and  i^t]  are  respectively  the  carrier  frequency  and  the  complex  modulation 
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function  of  the  radar  waveform,  is  a  scaling  factor  accounting  for  the  power  level 
transmitted  and  g^  is  the  complex  directivity  of  the  transmit  phase  center.  The 
constant  is  given  by: 


(29) 


where  Pj^  is  the  peak  transmitted  power. 

In  this  work,  i/[t]  is  a  train  of  rectangular  pulses  of  duration  Tp,  with  provision 
for  frequency  hopping  from  pulse  to  pulse: 

N  -1 
P 

t^t]  -  ^  n 

k-0 

where  /i[t]  is  the  single  pulse  modulation  function,  is  the  angular  frequency  increment 
for  the  k*^h  pulse  and  rr  is  the  Pulse  Repetition  Interval  (PRI).  The  definition  of  the 
rectangle  function,  n[t],  is 


t-kr 


L  ""P 


•/i|t  -kTjl-  expljw^t] 


(30) 


n[t]  =1  for  t  e  [0,1] 

=  0  elsewhere  (31) 

If  the  waveform  does  not  have  frequency  hopping  ,  uj-^  is  simply  set  equal  to  zero  in  Eq, 
30. 
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In  this  work,  we  allow  for  two  single  pulse  modulation  functions.  The  lirst  one, 
describes  a  waveform  without  frequency  modulation.  The  second  one, 
describes  a  waveform  with  linear  frequency  modulation  for  pulse  compression.  The 
equations  for  these  single  pulse  modulation  functions  are: 


and 


%fmW - 1 


(32) 


""  exp[j7rbt^]  0  <  t  <  Tp  (33) 

where  b  is  the  linear  frequency  sweep  rate  related  to  the  frequency  sweep,  B,  by: 


b  =  ±B/Tp  (34) 

The  sign  ±  in  front  of  B  indicates  that  the  frequency  of  the  RF  pulse  can  either  be 
up— shifted  or  down— shifted,  according  to  the  sign  chosen.  Since  B  is  approximately 
equal  to  the  inverse  of  the  compressed  radar  pulse  width,  r^,  we  have  that 

Ibl-l/rpT^  (35) 


2.4  The  Signal  from  a  Point  Source 

For  clarity  and  completeness,  we  describe  the  signal  at  various  points  along  the 
radar  receiver  shown  in  Fig.  3.  Let  yj[t]  be  the  signal  at  the  location  of  the  ith  phase 
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center.  This  signal  is  a  delayed  replica  of  the  transmitted  signal  backscattered  by  a 
point  target,  so  that: 


y[l‘)  =  g?  /)  u  I  -  Aj[t) 


(36) 


where  gj  is  the  directivity  of  the  ith  phase  center,  and  p  is  the  complex  reflectivity  of 
the  point  target.  The  constant  K  accounts  for  the  system  losses,  L  ,  and  the  geometric 
propagation  losses  and  is  given  by: 


K  = 
r 


A2 


(4x)2R4LsJ 


1/2 


(37) 


SAMP  ANALOis  TO  DIGITAL  CONVEPTFP 


Fig.  3  -  Radar  system  block  diagram 
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Normally  the  receiver  does  not  operate  when  the  transmitter  transmits  the 
waveform.  Therefore,  in  order  to  obtain  y^[t],  the  signal  actually  sensed  by  the  ith 
phase  center,  one  has  to  multiply  yj[t]  by  an  eclipsing  function,  E[t]; 

yfM  =  y'lt)  •  E[tl  (38) 

The  form  of  E[t]  can  be  determined  by  inspection,  from  Fig.  4. 


T 

r 


Fig.  4  —  Diagram  showing  the  transmitted  and  the  received  waveforms:  4a  shows 
the  transmitted  waveform  and  the  the  receiver  switching  time;  4b  is  the  eclipsing 
function;  4c  is  the  waveform  without  eclipsing;  4d  and  4e  show  eclipsing  of  respectively 
the  end  and  the  front  of  the  radar  pulses. 
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Fig.  4a  shows  the  transmitted  waveform.  We  assume  that  the  radar  transmits  the 
RF  pulses  on  a  continuous  basis  rather  than  transmitting  Np  pulses  and  turning  off  for 
some  time  before  repeating  the  process  over  again.  As  shown  in  Fig.  4a,  the  receiver  is 
inoperative  for  a  time  equal  to  t  12.  This  is  the  time  during  which  the  system  switches 

o 

from  the  transmit  to  the  receive  state.  The  time  during  which  the  receiver  is  functional 
is  shown  in  Fig.  4b.  This  is  the  eclipsing  function,  E[t].  It  is  a  series  of  rectangle 
functions  of  duration  given  by 


E[t|  =  X 

k=0 


where  is 


(39) 


(40) 


Fig.  4c  shows  the  received  pulse  train  for  the  case  where  there  is  no  eclipsing.  Fig. 
4d  and  4e  show  the  received  pulse  train  for  the  cases  where  the  front  and  the  end  of  the 
pulses  are  eclipsed.  Thus,  the  eclipsed  signal  is  a  train  of  pulses  with  each  pulse  being 
equal  to  or  shorter  than  Let  n  be  the  duration  of  the  eclipsed  pulses  and  let  be 
the  portion  of  the  leading  edge  of  the  pulse  which  may  have  been  eclipsed.  We  can  thus 
write  y^[t]  as: 
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yf[t]  =  K  p  exp  j£j^(t-Aj[t]) 


N  -1 
P 


I" 

k=0 


t-krj.-Aj[t]  exp|^jt^(t-Ajt])  (41) 


where  K=K,K  . 

V  1 

Let  k'  be  the  number  of  PRI’s  elapsed  before  the  first  pulse  from  the  target 
reaches  the  receiver.  The  parameter  k'  is  thus  largest  integer  such  that 

k'  <  (^)|T^  (42) 

To  calculate  the  values  of  and  Tj,  we  must  consider  four  possibilities. 

1)  No  eclipsing. 

This  occurs  when 


k'Tj+T  +Tj/2  <  2R/c  <  (k'  +  l)T^-T  -Tj/2 


(43) 


and  in  this  case 


Tj  =  0  (44) 

rj  =  Tp  (45) 
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2)  The  leading  edge  of  the  pulse  is  edipsed. 
This  occurs  when 


2R/c  <  k'r  +T  +T  /2 
r  p  s' 


and  in  this  case 


T,  =  T  —  Tf 

1  p  f 

3)  The  trailing  edge  of  the  RF  pulse  is  eclipsed. 

This  occurs  when 


2R/c  >  (k'+l)rj.-rp-Tg/2 


(46) 


(47) 

(48) 


(49) 


and  in  this  case 


Tf  =  0  (50) 

Tj  =  (k'  +  l)rj.-2R/c-rg/2  (51) 


4)  The  whole  RF  pulse  is  eclipsed. 

This  occurs  when 

k '  <  2R/c  <  k '  rj,+  tJ2  (52) 

and  in  this  case,  no  energy  is  detected  by  the  radar. 
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Once  the  signal  is  received  from  the  target,  it  is  amplified  and  converted  to 


baseband.  The  video  signal  is  thus: 


yjW  =  6^“’’  y?  W  exp 


-j“c‘ 


(53) 


where  is  the  gain  in  the  amplification  chain  including  Sensitivity  Time  Control 

(STC). 


To  a  very  good  approximation,  one  can  assume  that  for  the  envelope  of  the 
received  signal,  A.[t]  is  equal  to  the  gross  time  delay,  2R/c,  so  that,  using  Eq.  41  for 
y^[t],  the  final  expression  for  the  video  signal  is 


gj  P  exp 


X 


N  -1 
P 


I" 

k=0 


t— 2R/c-kr^— Tj 


t-krj.-Aj(t]  exp  jwj^(t-Aj[t]) 


(54) 


where 

Si  =  S  K,  st  5™’’  (55) 

In  the  next  sub  section,  Eq.  54  is  evaluated  separately  for  the  NFM  and  the  LFM 
waveforms. 
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2.4.1  The  NFM  Waveform 


As  shown  in  Fig.  3,  the  first  operation  performed  on  the  NFM  signal  is  match 
filtering.  In  general,  the  impulse  response  of  a  filter  matched  to  a  waveform  f[t]  is  given 
by: 


hit)  =  K„,  (56) 

where  is  a  normalization  constant.  In  this  work,  the  received  signal  is  matched  to  a 
single  rectangular  pulse  of  the  transmitted  burst.  However,  the  mathematical 
expression  for  the  matched  filter  may  or  may  not  include  the  frequency  steps  To 
account  for  these  two  possibilities,  the  center  frequency  of  the  matched  filter  is  shifted 
by  where  wj^=0  for  clutter  suppression  and  for  ECCM. 

To  simulate  the  operation  of  a  clutter— lock  circuit  (Ref.  6),  the  transfer  function  of 
the  matched  filter  is  translated  by  the  doppler  frequency  o;^.  In  this  manner,  its  center 
frequency  corresponds  to  the  Doppler  shifted  frequency  of  main  beam  clutter.  The 
impulse  response  of  the  matched  filter  for  the  NFM  waveform  is  thus: 


where 


^m^^l  “  ^NFM  ^1  ^/’’p 


exp 


j"dk‘ 


(57) 


“dk  -  “d  +  "k 


(58) 


and  where  for  the  NFM  waveform.  Equation  (57)  implies  that  the  width 
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of  the  range  gates  is  equal  to  Tp.  In  practical  systems,  this  is  not  always  the  case.  To 
account  for  this  possibility,  the  impulse  response  of  the  NFM  filter  is  written  as 


''NFmI'I  "  ’^NFM 


(59) 


where  r  is  the  width  of  the  range  gates  (r  <  r  <  r  ). 
g  o' 


Requiring  that  the  transfer  function  has  unit  amplitude  at  w  =  we  find  that 
^NFM 


^NPM  ~  ^/^g 


(60) 


The  matched  filtered  signal,  z.(t],  is  obtained  by  convolving  Eq.  54  with  Eq.  59  to 


give: 


Zi[t]  =KNFMSi^  I  expl-jWpAjN) 


— OD 


N  -1 
p 


X" 

k-0 


Tf—2Rf  c-kr^—T^ 


^1 


exp 


n 


_(^JexpJjw^j^(t-77)j  dr? 

s 


(61) 


Defining 


T}'  =  rj-t 


(62) 
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and 


a[t]  =  — t+2R/c+krj,+r£ 


(63) 


we  obtain 


z-[t]  —  p  exp 


X 


N  -I 

P  QD 

^  rexpJ-j(a;^+a;j^)A.[7/'+t]  +  « 

k=0 

—00 


n 

1 

N 

1 

_ ! 

n 

L  J 

■V 

(64) 


As  shown  in  Fig.  5,  the  product  of  the  two  rectangle  functions  in  Eq.  64  differs 
from  zero  only  over  an  interval  defined  by 


(65) 


Inequality  65  can  be  rewritten  as: 


t— 2R/c— Tp  <  kr^  <  t— 2R/c— r^+Tg 


(66) 


Since  T^-¥r^<T^,  Eq.  66  shows  that  for  any  given  values  of  t  and  R,  there  is  only  one 
value  of  k,  Ic,  satisfying  the  inequality.  Physically,  this  means  that  only  the  pulse  in 
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In  equation  67  and  68,  the  value  of  r  is 


=  TgH  a[t]  +rjH|^— a[t]j  for  a[t]<0  or 

a[t]>0  A  rj+Q![t]>rp 

r=  T^+a[t]  otherwise 


Where  H[t]  is  the  Heaviside  step  function  given  by: 


H[t]  =  1  for  t  >  0 

=  0  elsewhere 


Hence,  making  use  of  Eq.  67  —  70  and  replacing  A.[t]  by  its  value  Eq.  24—27,  Eq 


64  becomes 


Zi[t]=KNFM^^  exp  jw.t  K 

1' 

max 

J  exp  ^Q+^j(7/'+t)+y^'+t)^  +j^dk^' 


if  a  value  of  k  in  the  interval  [0,Np— 1]  can  be  found.  Otherwise,  Eq.  64  is  equal  to  zero 
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In  Eq.  71,  and  are  respectively 


n 


ck  =  +  ‘^c 


and 


n 


dk  =  “dk 


(72) 

(73) 


Since  the  values  of  and  do  not  exceed  10“^s  for  typical  ABRs  and  since 
2 

^2  is  small,  the  terms  62^'  and  62^'^  are  negligible.  Extracting  the  terms  independent 
of  T}'  from  the  integral  we  obtain  after  performing  the  integration 


z.[t]  = 


2gi/7 

"nr 

k 


r  f  _  .  .  .11 

-jnc  (7-+a[t])' 

1 

0 

exp 

iw 

2 

sin 


where 


^k  “  ^ck^l  “  ^dk 


(74) 


(75) 


The  tilda  symbol  in  Eq.  74  is  reminder  that  the  corresponding  parameters  depend 
on  a  unique  value  of  k=k. 

To  simplify  the  notation  for  the  later  calculations  of  the  clutter  covariance  matrix 
in  Sect.  2.6,  we  extract  all  the  terms  independent  of  p  in  the  expression  of  z.[t]: 


Zj[tl  =  p  z.[t] 


(76) 
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where 


exp 


(r+a[t]) 


sin 


2 


(77) 


In  modern  radar  receivers,  an  analog  to  digital  converter  (ADC)  samples  the  signal 

for  input  to  the  digital  signal  processor.  Using  Eq.  74,  one  obtains  the  1*^  sample  of  the 

signal  for  the  range  gate  at  a  distance  R  from  the  radar  by  letting  t  equal 

6 


t  = 


(78) 


In  the  last  equation,  1  is  an  element  of  the  interval  [0,Np-l]. 
2.4.2  The  LFM  Waveform 


With  an  LFM  waveform,  it  is  no  longer  practical  to  implement  a  matched  filter  to 
process  the  signal.  This  is  because  the  transfer  function  of  the  matched  filter  contains 
ripples  which  are  difficult  to  implement  in  practice.  This  applies  especially  when  the 
compression  ratio  is  smaller  than  about  twenty  (Ref.  7,8).  Moreover,  the  compressed 
pulses  exhibit  large  sidelobes  which  must  be  weighted  down  to  reduce  the  range 
sidelobes.  It  is  thus  common  practice  to  implement  a  mismatched  filtered  (MMF)  with 
LFM  waveforms.  We  thus  follow  Ref.  7  and  model  a  MMF  consisting  of  compressing 
the  received  signal  and  weighting  it  in  the  frequency  domain.  The  transfer  function  of 
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the  compressor  is 


H^[w]  =  exp 


j(“  -  “dk)'' 

47rb 


(79) 


As  for  the  NFM  case,  the  center  frequency  of  the  filter  is  shifted  by 


One  can  easily  verify  that  H^[a;]  is  a  compressor  by  using  the  group  delay  concept. 
The  group  delay  of  a  filter,  r^,  relates  to  the  velocity  at  which  information  or  energy 
flows  in  a  system.  It  is  given  by  (Ref.  8) 


(80) 


where  0[u)]  is  the  phase  of  the  transfer  function  of  the  filter.  Using  Eq.  79,  the  group 
delay  for  the  compressor  is; 


("-"dk) 

^d  27rb 

Putting  Wjjj=0  to  simplify  the  interpretation,  this  equation  shows  that  the  group  delay 
is  proportional  to  the  angular  frequency  of  the  waveform.  Hence,  on  transmit,  the 
frequency  of  the  waveform  is  linearly  increased  at  a  rate  b.  On  receive,  the  compressor 
decreases  it  at  the  same  rate. 

For  the  sake  of  simplicity,  we  use  an  exponential  weighting  of  the  frequency 
spectrum  of  the  waveform  to  achieve  range  sidelobe  reduction  (Refs.  9,10).  Although  it 
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might  not  be  the  preferred  weighting  function  from  a  practical  point  of  view,  it  performs 
like  the  more  usual  weighting  schemes.  The  weighting  operation  can  be  represented  by 


Hp(a;)  =  exp 


1  -  “dfc)'' 

47r 


(82) 


where  q  is  the  weighting  coefficient.  It  is  adjusted  to  achieve  an  attenuation  of  A^g 
decibels  at  the  edges  of  a  frequency  window  extending  B  Hz  on  either  side  of  the  signal 
center  frequency.  One  can  easily  show  that  q  is: 


q  =  0.147  Ajg/B^  (83) 

The  overall  impulse  response  of  the  MMF  is  obtained  by  computing  the  Fourier 
transform  of  the  product  of  Eq.  79  and  Eq.  82.  Using  the  work  of  Prudnikov  et  al.  (Ref. 
11)  the  result  is 


where 


‘LFM 


[t]  =  K 


LFM 


exp 


t^ 


q"b"  +  l 


(84) 


^LFM  “ 


44^ 

jb  q^  +  1 


(85) 


Eq.  84  shows  that  the  impulse  response  of  the  MMF  is  centered  on  t=0  rather  than 
on  r  /2  for  the  NFM  case.  Hence,  for  the  LFM  calculations  alone,  it  is  convenient  to 

O 
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redefine  the  rectangle  function  in  the  following  manner: 

n  [t]  =  1  for -1/2  <  t  <  1/2 

=  0  elsewhere  (86) 

In  this  manner,  the  maximum  output  of  the  filter  occurs  when  the  range  to  the  target, 

R,  equals  R  ,  as  is  the  case  for  the  NFM  waveform. 

S 

Following  the  procedure  used  for  the  NFM,  the  signal  at  the  output  of  the  MMF  is 
the  convolution  of  yjt]  with  h[t]: 


N  -1 
P 


n 


k=0 


z.[t]  -  g.pj  exp 

— O) 

Tf-2R/ c-kr^—T^ 


exp 


j7rb 


exp 


exp 


d?; 


(87) 


Defining 


and 


T]'  —  T]  —  2R/c  —  kXj,  — 


(88) 

(89) 
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and  using  Eq.  63,  Eq.  87  becomes 


V‘l  =’^LFMSi'’“P 


N  -1 
P 


r,/2 


exp 


exp  [h4]  / 

-rj/2 

2I _ [  xb^j  +  bq)*' 


k=0 


j7rb  7/'+2R/c+rj>-Aj 


exp 


q^b^+1 


Ti'  +  a[t] 


drj' 


(90) 


In  this  equation,  the  complex  exponentials  are  weighted  by  a  series  of  Gaussian 
functions  (the  second  term  of  the  second  exponential  function): 


exp 


xqb"^ 

qV+1 


r/'-t+2R/c+kTj.+r£ 


12 


(91) 


The  first  of  these  Gaussian  functions  (k=0)  is  centered  at 


77'  =  t— 2R/c— Tj. 


(92) 


and  the  others  are  r^,  apart  around  t}'  .  The  width  of  the  Gaussian  functions  between  the 
point  Tj'=0  and  the  e~^  point,  K,  is  by 


2.2^; 

q  b  -H 


xqb ' 


1/2 


(93) 
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The  parameter  H  can  be  expressed  in  terms  of  r  and  r  for  more  clarity.  Since  q  is 

c  p 

2 

proportional  to  (Eq.  83)  and 


we  find  that 


T_T 


p  C 


H  oc  r 

P 


(94) 


(96) 


since  tp 


is  normally  at  least  ten  times  greater  than  t^. 


It  is  interesting  to  note  that  the  width  of  bj^pj^[t]  is  proportional  to  although 
the  width  of  hp[t]  is  proportional  to  r^.  The  width  of  hppj^[t]  has  to  be  proportional  to 
Tp.  Otherwise,  only  a  small  fraction  of  the  energy  returned  by  the  target  would  be 
processed  by  the  receiver.  This  would  result  in  a  low  the  signal  to  noise  ratio. 


The  interval  of  integration  in  Eq.  90  extends  only  over  [-rj/2,rj/2].  It  is  thus 
likely  that  one  of  the  Gaussian  functions  contributes  significantly  more  than  the  others 
to  the  integral.  This  can  be  confirmed  with  the  help  of  the  Error  function  (Ref,  12). 
The  first  case  of  interest  is  when  one  of  the  Gaussian  overlaps  the  interval  of 
integration.  Assuming  a  worst  case  duty  cycle  of  20%  and  no  eclipsing,  the  power 
contribution  of  the  Gaussian  centered  on  the  interval  of  integration  is  more  than  100  dB 
larger  than  the  contribution  from  the  adjacent  Gaussians.  The  second  case  of  interest  is 
when  the  Gaussian  functions  are  symmetrically  located  on  either  side  of  the  interval  of 
integration.  Their  power  contribution  to  the  integral  is  approximately  47  dB  lower  than 
that  of  the  Gaussians  centered  on  the  interval  of  integration.  These  numbers  show  that. 
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like  for  the  NFM  case,  it  is  sufficient  to  consider  the  contribution  of  only  one  Gaussian 
to  the  integral  of  Eq.  90.  This  is  the  Gaussian  which  is  located  the  closest  to  the  center 
of  the  integration  interval.  The  index  of  this  Gaussian,  k,  is  the  integer  element  of 
[0,Np— 1]  which  minimizes  the  expression  t— 2R/c— kr^— r^. 

Expanding  the  value  of  Aj[t]  and  assuming  that  only  the  k^^  pulse  contributes  to 
the  integral,  Eq,  90  becomes 


ZiN  =  KfiEM  Si  ^  "j^^dk^W  exp  x 

rj/2 

-rj/2 

exp  j7rb|7?'+2R/c+rp(5Q-<5j  -^2  ^'+^k  ^ 

r  r  1 21 

exp  -  r?'  +  a[t]  dr;'  (96) 

q  b^+l[  J 

Using  Table  I,  one  finds  that  many  terms  in  Eq.  96  are  negligible  over  the  interval 
of  integration.  The  last  equation  thus  reduces  to 

Zi[t]  =  gj  pexp  -ji^(^j^a[t]  exp 

exp [/Jji] -  ^ 

L  L  J  q  Ij  ^.1  u  j  j 
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/  “p[-j{“ck*r”dk  +j 

-r,/2 


exp  jjjb  -  ’T?'^  dp' 


•2  -  _;,K  i  xMi+bi 


X'^  =  — j^rb  + 


q  b  +1 


Y=-(l/2X^)|jfi^j^<5^  -  -  j27rb 


Trb(. i  +  bq^ 
q^b^+1 


By  completing  the  squares  we  obtain  after  some  manipulations 


Zi[t]  =  -^^-g-pexp  -j^^^i^aft]  exp 

. 


exp  -j^^ck^i  -  -  -^^y^a^[t]  +  jTrb  ^  * 

L  L  J  q  b  +1  *■  J  J 


,/?/2|Erfc  X  Y-T|/2  -  Etfc  X  Y+T|/2 


(100) 


where  Erf(x]  is  the  complementary  Error  Function  defined  as 


Erfc[x]  =  1  —  Erffx] 


(101) 
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In  order  to  simplify  the  notation  for  the  calculation  of  the  clutter  covariance 
matrix  in  Sect.  2.6,  the  terms  independent  of  p  are  extracted  out  of  the  expression  for 
z.[t]: 


z.[t]  =  p  z.[t] 


(102) 


where 


exp 


exp 


-j^ck^i 


Pi 


_  Ml  +  .bq)a^[t]+jxb[r^/?r,^  ] 
q^b^  +  1  L  1  ^  ij 


Vi/2|Erfc  X  Y-r|/2jl  -  Erfc  X  Y+T|/2 


(103) 


Eqs.  76  and  102  describe  the  signal  received  from  a  point  target  for  respectively 
the  NFM  and  LFM  waveforms.  They  can  also  be  used  to  obtain  the  signal  from  a 
distributed  target.  Our  interest  in  this  work  is  primarily  ground  clutter  against  which 
point  targets  have  to  compete  for  detection.  This  is  the  subject  of  the  next  two  sections. 
The  following  calculations  can  easily  be  extended  to  volumetric  clutter  (rain  or  chaff). 
A  similar  procedure  would  also  apply  to  compute  the  signal  from  complex  targets  that 
can  be  resolved  by  the  radar. 


2.5  The  Signal  from  a  Distributed  Target 
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The  signal  from  a  distributed  target  is  the  sum  of  the  signals  from  a  collection  of 
contiguous  elemental  scatterers  spanning  the  target.  In  the  case  of  the  clutter,  the 
reflectivity  of  an  elemental  clutter  patch  of  area  dA^  is 

Pc  = 

In  Eq.l04,  p°[7j  is  the  clutter  reflectivity  per  unit  area  and  is  a  function  of  the  grazing 
angle,  7^.  Ref.  1  shows  that  the  area  of  a  clutter  patch  on  a  spherical  surface  (the  earth) 
is 

where  the  subscript  c  refers  to  a  clutter  patch.  The  signal  from  an  elemental  clutter 
patch,  dcjj,  is  given  by  Eqs.  76  and  77  for  the  NFM  case  or  by  Eqs.  102  and  103  for  the 
LFM  case.  For  clutter,  replaces  p  and  the  range  and  the  velocity  parameters  are 
those  of  the  clutter  element  under  consideration. 

The  total  clutter  return  is  the  integral  of  the  elemental  contributions  over  the 
terrain  within  the  footprint  of  the  antenna; 

c,(tl  =  /  /  dCjUl  (106) 
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The  limits  of  integration  extend  at  most  over  the  following  interval; 


and 


t  [0,21) 


(107) 


fK,(r!-rV/2 


(108) 


2  2  1/2 

where  h  is  the  altitude  of  the  aircraft  and  (r„  —  r  )  '  is  the  range  to  a  point  on  the 
Earth’s  horizon  as  seen  from  the  radar.  Clearly,  because  of  the  weighting  effect  of  the 
radar  antenna  gain  pattern,  only  a  small  fraction  of  the  clutter  patches  contribute 
significantly  to  the  clutter  return. 


2.6  The  Clutter  Covariance  Matrix 


As  discussed  in  Refs.  1^,  the  clutter  covariance  matrix,  C,  is  required  to  compute 
the  average  interference  power  against  which  the  target  is  competing  for  detection.  In 
the  general  case,  C  is  a  measure  of  the  correlation  between  the  clutter  signals  found  on 
subapertures  i  and  i'.  The  components  of  C,  c-^,,  1  e  [0,Np— 1]  are 


(109) 


where  <>  is  the  ensemble  average  operator.  In  Eqs.  76  and  102,  the  only  statistical 
parameter  is  p^.  As  discussed  in  Refs.  1-4,  the  statistics  of  p^  are  appropriately 
described  by 


(110) 
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and 


where  (5{x— x']  is  the  Dirac  delta  function.  The  first  expression  states  that  the  mean 
value  of  the  complex  reflectivity  of  a  clutter  patch  is  zero.  The  second  expression  states 
that  distinct  clutter  patches  are  uncorrelated  and  that  the  backscattering  coefficient  of  a 
clutter  patch  is 

Using  Eqs.  110  and  111  and  noting  that  the  statistical  independence  of  the  clutter 
patches  causes  the  cross  terms  to  vanish,  Eq.  109  becomes: 

Eq.  112  can  easily  be  expanded  by  replacing  Cj[t]  and  C|,[t]  by  their  respective 
values.  The  expression  obtained  is  rather  complicavcJ  and  not  very  informative.  If 
frequency  hopping  is  not  used,  the  expanded  version  of  Eq.  112  simplifies  considerably. 
For  the  NFM  case,  Eq.  112  is  identical  to  the  equation  found  in  Refs.  1-4,  provided  the 
second  order  phase  terms  are  neglected. 


2.7  The  Thermal  Noise  Covariance  Matrix 


In  radar  systems,  the  target  also  has  to  compete  against  white  thermal  noise  for 
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detection.  As  for  clutter,  the  noise  covariance  matrix,  N,  is  used  to  compute  the  noise 
interference.  Its  definition  is 


N(l,r]  =  ^n|ljn  (l+r]^ 


(113) 


where  r  is  the  time  interval  and  where  n[t]  denotes  the  noise  signal.  Assuming  that  the 
noise  is  a  stationary  stochastic  process,  N  becomes  independent  of  t: 


N(r|  =  <J^n(t|n*(t+rl^ 


(114) 


As  discussed  in  Ref.  13,  the  expected  value  in  the  expression  above  is  giver  by  the 
Weiner— Kinchine  relation: 

QD 

N[r]  =  J S^[fl  exp[j27rfr]  df  (115) 

— GD 


where  f  is  the  frequency  and  Sj^[f]  is  the  autospectral  density  function  of  the  noise.  At 
the  output  of  a  filter  of  transfer  function  H[f),  S^[{]  relates  to  the  autospectral  density 
function  of  the  white  noise  at  the  input  of  the  filter,  S  [f],  in  the  following  manner: 

=  (“«) 

For  white  thermal  noise,  S^[f]  is  a  constant  and  Nq  is  given  by  the  well  known 
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relationship 


No  =  keTs 


(117) 


where  kg  is  the  Boltzman  constant  and  is  the  equivalent  system  temperature.  Using 
Eq.  59  to  compute  replacing  into  Eq.  116  and  using  Ref.  11,  we  find 

N[r]  =  Nq  (Tg  -  r)  for  r  < 

N[r]=  0  elsewhere  (118) 


Replacing  r  by  (1-1 ')rj.  and  by  ‘Is  value  (Eq.  60),  we  obtain 


(119) 


where  6\\'  is  the  Kroenecker  delta  function.  Eq.  119  is  the  well  known  relation  linking 
the  noise  power  to  the  system  temperature  and  the  bandwidth  of  the  radar  signal 
processor.  Finally,  assuming  that  the  noise  is  uncorrelated  from  one  receive  phase 
center  to  another,  Eq.  119  becomes 


(120) 


Using  Eqs.  84  and  85,  the  equivalent  expression  for  the  LFM  signal  is 


kg  Tj/A 


‘*11" 


^.i' 


(121) 


Again,  since  qotr 


2 

c’ 


the  noise  power  is  proportional  to  the  bandwidth  of  the  radar  signal. 
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2.8  The  Signal-to-Interference  Ratio 

The  point  target  competes  against  both  clutter  and  noise  for  detection.  The  total 
interference  covariance  matrix,  A,  describes  the  statistical  properties  of  the  combined 
interference.  Since  the  clutter  and  the  noise  are  statistically  independent,  A  is 

A  =  C  +  N  (122) 

where  C  and  N  are  respectively  the  clutter  and  noise  covariance  matrices. 

Recalling  the  block  diagram  of  Fig.  2,  after  analog  to  digital  conversion  is 
performed,  a  digital  signal  processor  is  used  to  enhance  the  target  signal  against  the 
interference.  It  also  extracts  some  features  of  the  target.  This  enables  the  data 
processor  to  provide  useful  information  about  the  target,  such  as  position  in  space, 
relative  velocity  etc.  The  action  of  the  signal  processor  is,  in  the  ideal  case,  a  linear 
filtering  operation.  Mathematically,  the  filtering  operation  is  the  scalar  product  of  the 
signal  vector  by  a  weight  vector,  W.  The  components  of  W  are  w^j,  where  i  is  the 
receive  phase  center  index  and  1  is  the  sample  number.  At  the  output  of  the  signal 
processor,  the  ratio  of  the  target  power  to  the  interference  power  is  the 
Signal-to-Interference  Ratio  (SIR): 

SIR  =  JsLlW__L  (123) 

W '  •  A  •  W 

where  *  is  the  complex  conjugate  operation,  T  is  the  transpose  operation  and  f  is  the 
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complex  conjugate  and  transpose  operations.  Since  S  is  a  function  of  the  velocity  of  the 
target,  SIR  also  varies  with  v^. 

The  denominator  of  Eq.  123  contains  the  average  interference  covariance  matrix. 
Hence,  SIR  is  an  average  value.  The  approach  adopted  in  this  report  allows  one  to 
obtain  the  average  performance  of  the  radar  system  without  using  Monte  Carlo 
techniques. 


3.0  NUMERICAL  RESULTS 


3.1  The  Computer  Simulation 

The  computer  simulation  developed  for  SBRs  in  Ref.  4  was  modified  to  apply  for 
the  airborne  scenario.  The  main  modifications  consisted  of  implementing  the  equations 
developed  in  this  work,  decoupling  the  velocity  of  the  radar  platform  from  the  satellite 
orbital  dynamics,  allowing  the  aircraft  to  dive,  implementing  an  air— to— air  mode  and 
removing  the  Earth’s  rotation  since  the  aircraft  navigates  with  respect  to  the  surface  of 
the  Earth.  An  upcoming  report  will  describe  the  computer  simulation  in  detail. 

3.2  The  Baseline  Airborne  Radar  System 

Table  I  describes  the  baseline  Airborne  Intercept  Radar  system.  It  is  a 
conventional  pulse  Doppler  radar  using  a  Medium  PRF  (MPRF)  waveform  (Ref.  14). 
For  the  baseline  case,  the  aircraft  flies  at  500  m/s  at  an  altitude  of  5  km. 
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TABLE  m 


AIRBORNE  INTERCEPT  RADAR  PARAMETERS 


Power 

1.0  kW 

Pulse  Repetition  Frequency 

9.1  kHz 

Compression  Factor 

13.0 

Range  Resolution 

250  m 

Number  of  Pulses  Transmitted 

16 

Antenna  Type 

Circular 

Antenna  Diameter 

0.7  m 

Aperture  Illumination  Function 

Taylor 

Antenna  Sidelobe  Level 

^0  dB 

Platform  Velocity 

500  m/s 

Platform  Altitude 

5  km 

Target  Altitude 

0  m 

Cross  Section  Fluctuations 

SW  II 

Terrain  Type 

Land 

Signal  Processing 

DFT 

Time  Domain  Weighting 

Blackman 

In  the  first  example,  the  radar  is  in  the  down— looking  mode,  attempting  to  detect 
a  slow— moving  target  at  a  range  of  10  km.  The  target  is  close  to  the  surface  of  the 
Earth.  The  radar  transmits  an  LFM  waveform  with  a  compression  ratio  of  13.  After 
pulse  compression,  the  signal  is  coherently  integrated  using  the  standard  FFT  with 
Blackman-Harris  weighting.  Furthermore,  two  bursts  of  pulses  are  integrated 
non— coherently  to  reduce  the  effect  of  target  scintillation.  We  assume  that  the  signal 
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from  the  target  does  not  ^uctuate  during  coherent  integration.  The  target  cross  section 
fluctuations  obey  the  Swerling  II  model  from  coherent  burst  to  coherent  burst. 

Fig.  6  shows  the  signal-to-interference  ratio  (SIR)  of  the  radar  as  a  function  of 
the  velocity  of  the  target.  The  dotted  line  is  the  detection  threshold  above  which  a 
target  is  present.  Fig.  6  shows  that  the  SIR  of  the  radar  is  very  small  when  the  target  is 
stationary.  At  low  velocity,  the  signal  processor  cannot  distinguish  the  target  from  the 
clutter  on  the  basis  of  its  Doppler  shift.  As  it  goes  faster,  the  Doppler  shifted  target 
signal  moves  away  from  the  spectrum  of  the  mainbeam  clutter.  Eventually,  the  SIR 
exceeds  the  detection  threshold  and  the  target  is  detectable.  In  the  flat  portion  of  the 
SIR  curve,  the  clutter  is  almost  completely  filtered  out  and  the  radar  performance  is 
essentially  noise  limited.  In  the  regions  of  very  low  SIR,  the  radar  performance  is 
essentially  clutter  limited.  Because  sampled  signals  are  periodic  with  a  period  equal  to 
the  PRF,  the  SIR  becomes  very  small  again  when  the  target  moves  at  approximately 
150  m/s. 

Hence  for  the  baseline  system,  the  target  is  detectable  in  the  following  velocity 
intervals  ;  (-136,-18  m/s)  and  (15.5,133.5  m/s).  This  means,  for  instance,  that  the 
radar  is  unlikely  to  detect  ground  vehicles  unless  they  have  a  radial  velocity  exceeding 
56  km/h.  This  corresponds  to  a  ground  velocity  of  65  km/h,  assuming  an  angle  of  30° 
between  the  velocity  of  the  aircraft  and  the  antenna  ajds  of  symmetry.  The  asymmetry 
between  the  positive  and  the  negative  detectable  velocity  interval  is  attributed  to  the 
range-doppler  coupling  of  the  LFM  waveform. 
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Fig.  6  —  The  Signal— to-Interference  Ratio  for  the  baseline  system 

There  is  no  trace  of  an  altitude  return  (the  return  from  terrain  directly  under  the 
aircraft)  in  Fig.  6.  The  altitude  of  the  radar  being  5  km  and  the  ambiguous  range  being 
16.48  km,  the  altitude  return  occurs  beyond  a  range  of  21.48  km.  Fig.  7  compares  the 
SIR  curves  of  the  baseline  case  with  the  one  for  a  range  equal  to  21.48  km.  The  altitude 
return  causes  the  Detectable  velocity  intervals  to  practically  disappear.  Hence,  the 
radar  could  hardly  detect  a  10m2  target  at  this  range. 
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Fig.  7  —  SIR  of  the  baseline  system  with  altitude  return 

Fig.  6  and  7  show  that  the  Intervals  of  Detectable  Velocity  shrink  as  the  PRF 
decreases  since  the  areas  of  small  SIR  move  closer  to  each  other.  A  radar  with  a  PRF 
equal  to  or  lower  than  6  kHz  would  be  ineffective  in  this  situation.  If  the  PRF  of  the 
radar  is  too  large  on  the  other  hand,  the  range  ambiguities  increase  and  the  target  has 
to  compete  with  more  ambiguous  clutter  for  detection.  Thus,  a  MPRF  waveform  for 
which  the  radar  signal  is  ambiguous  in  both  range  and  Doppler  is  a  good  compromise  for 
the  detection  of  moving  targets  in  the  look— down  mode. 
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The  performance  of  the  radar  with  an  LFM  waveform  can  be  approximated  quite 
closely  by  using  an  NFM  waveform.  In  this  case,  the  pulse  width  is  scaled  down  by  the 
compression  ratio  and  the  peak  power  is  scaled  up  by  the  same  amount.  This  is  shown 
in  Fig.  8  where  the  SIR  for  the  LFM  waveform  is  compared  to  the  "equivalent"  NFM 
waveform.  In  the  noise  limited  region,  the  SNR  is  approximately  3  dB  larger  for  the 
NFM  case  than  for  the  LFM  case.  This  is  because  the  NFM  waveform  does  not  require 
weighting  for  range  sidelobes  reduction.  The  NFM  curve  is  also  slightly  narrower 
because  it  does  not  have  range  sidelobes.  Since  the  NFM  waveform  requires  much  less 
computer  time,  it  is  often  useful  to  make  this  approximation. 
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Fig.  8  —  SIR  of  the  LFM  and  the  NFM  "equivalent"  waveform 
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In  an  air— to— air  encounter,  the  main  beam  of  the  antenna  may  or  may  not  be 
pointed  toward  the  surface  of  the  Earth.  If  it  is,  Earth  clutter  penetrates  the  system 
and  interferes  with  target  detection.  However  if  the  main  beam  of  the  antenna  does  not 
intersect  the  surface  of  the  Earth,  clutter  penetrates  the  radar  system  only  through  the 
sidelobes  of  the  antenna.  Except  for  range  gates  in  which  the  altitude  return  is 
significant,  the  clutter  power  is  generally  quite  small  and  the  performance  of  the  system 
is  essentially  noise  limited.  This  is  shown  in  Fig.  9  in  which  the  baseline  radar  system  is 
looking  at  a  target  flying  at  an  altitude  of  7  km.  There  is  very  little  clutter  interfering 
with  the  target.  The  system  is  essentially  noise  limited. 
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Fig.  9  -  SIR  of  the  baseline  system  for  an  air-to-air  encounter 
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4.0  CONCLUSION 


In  this  report,  we  have  developed  a  theory  to  model  the  operation  of  airborne 
radars.  It  simulates  multiple  waveforms  and  multiple  phase  centre  systems  with 
frequency  hopping.  It  can  assess  the  performance  of  airborne  radar  systems  and 
investigate  advanced  interference  suppression  techniques.  A  very  interesting  feature  of 
the  model  is  its  ability  to  evaluate  the  average  performance  of  the  radar  without  the  use 
of  Monte  Carlo  techniques.  The  model  has  been  applied  to  a  baseline  Airborne  Intercept 
radar  to  demonstrate  its  capability  and  validity. 
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